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Abstract. We provide a new analytical approach to operator splitting for 
equations of the type ut = Au + uux where A is a linear differential operator 
such that the equation is well-posed. Particular examples include the viscous 
' Burgers' equation, the Korteweg-de Vries (KdV) equation, the Benney-Lin 

^ equation, and the Kawahara equation. We show that the Strang splitting 

, method converges with the expected rate if the initial data are sufficiently reg- 

' ular. In particular, for the KdV equation we obtain second-order convergence 

' in for initial data in ff'+s with arbitrary r > 1. 

1. Introduction 

I In this paper we study initial value problems of the type 

• ; (1.1) Ut = P{d.^)u + uu^ {xeR, 0<t< T), u\t^to = uo, 

, with a polynomial P of degree £ >2 satisfying 

^' (1.2) ReP{i()<0 for all ^gM. 

This class includes several important equations like 

Ut — Uxx + uux, the viscous Burgers equation, 

00 ' Ut — Uxxx + uux, the Korteweg-de Vries equation, 

^ . Ut ^ -Uxxx - P{uxx + Uxxxx) - Uxxxxx + uUx, the Benney-Lin equation [21 [H], 

. Ut — Uxxxxx — Uxxx + UUx, the Kawahara equation [12] . 

(N 
O ■ 

We employ operator splitting (see, e.g., [B]), i.e., construction of an approximate 
solution by concatenating the solutions of the separate problems 

(1.3) vt = P{dx)v and Wt = wwx- 
■ f-n 

^ \ More precisely, if we let u{t) = $^(uo) denote the solution of the initial value 

M ' problem 

(1.4) Ut = C(u), u\t=o = uo, 

then (sequential) operator splitting claims that the exact solution u{t) = <I>^(uo) is 
well approximated by w„, at t = nAt < T and as At — > 0, where 

(1.5) Un+i^^f{^%\u^)), n = 0,l,2,..., 
and C ^ A + B . In our case we have the operator 

(1.6) Au = P{dx)u 
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B{u) ~ UU; 



•X 



acting on appropriate Sobolev spaces. In the present paper a more refined operator 
splitting, known as Strang splitting, will be discussed. Strang splitting means that 
we consider the approximation 



The technique we use requires a well-posedness theory for the full equation in 
Sobolev spaces. This is currently available for the equations listed above, see [Tl] 
(the viscous Burgers equation), [U [T31 HH] (the KdV equation), jB^ (the Benney- 
Lin equation), and '111 (the Kawahara equation). The challenge in using operator 
splitting for the present class of equations is the fact that the Burgers step, i.e., 
solving the second equation in (|1.3p . generically introduces singularities in finite 
time irrespective of the smoothness of the initial data, while the solution of the 
full problem (jl.4p remains smooth. Thus the determination of the time step used 
in the operator splitting and the control of the smoothness of the approximate 
solution are rather delicate. This implies that the standard estimates that are 
required for operator splitting, see, e.g., [3], do not apply. The idea of applying 
operator splitting to the KdV equation was introduced in a short note by Tappert 
[19] ■ with numerical tests. A further study was done in [7] where a Lax-Wendroff 
theorem was proved, namely assuming convergence of operator splitting, one can 
prove that the limit is a weak solution of the KdV equation. Furthermore, extensive 
numerical computations were presented, all indicating that operator splitting does 
converge. A completely different and novel approach was taken in the paper [5]. 
Previous results had as goal to use operator splitting to show existence of solutions 
of the KdV equation. However, in [S] one took for granted the well-posedness of 
the Cauchy problem for the KdV equation in Sobolev spaces of high order. By 
using a bootstrap principle (see [18]), one could show that the Burgers step in a 
Sobolev space of lower order would avoid blow-up, thereby securing that the full 
approximation remained smooth, and indeed converged to the solution of the KdV 
equation. A key part of this analysis was to secure a uniform time step in the 
approximation. The approach in [5] also presented a new interpolation between the 
discrete time steps based on the introduction of a new auxiliary time variable, and 
hence a two- variable approximation was used. 

In the present paper we take a different approach. We do take as a starting point 
the well-posedness of the Cauchy problem for the KdV equation, and we use the 
same analysis to estimate carefully the change in Sobolev norm during the Burgers 
step. However, instead of introducing a two-variable approximation and applying 
the bootstrap principle, we perform an analysis that identifies the principal error 
terms of the local error as quadrature errors, which are estimated via bounds of 
Lie commutators. Such a type of analysis was first done for operator splitting 
for linear evolution equations in jlO) and for non-linear Schrodinger equations in 
[T7] . The error propagation is studied via "Lady Windermere's fan" [5], as is 
usual in the error analysis of time discretization methods, controlling here carefully 
the spatial regularity of the numerical approximations by first establishing lower- 
order convergence in a higher-order Sobolev norm, as in [T7]. In this way we 
improve and simplify the results in [5]. The main result says that the Strang 
splitting converges to the full solution with the expected second-order convergence 
rate, in the Sobolev norm for r > 1 in the case of H^+^ initial data for the 
KdV equation. Furthermore, as was clear already in [5], the analysis extends to a 
class of dispersive equations, and here we find it suitable to study the class (jl.ip 
subject to the condition (|1.2p . This means that the equations all share a Burgers 



(1.8) 



n = 0,1,2,... . 
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nonlinearity. The extension of the present approach to dispersive equations with 
other nonhnearities, not necessarily quadratic and also in higher spatial dimensions, 
appears feasible as long as an analogue of the regularity properties of the Burgers 
step can be established. We also refer to fS. where the approach of is applied to 
the generalized active scalar equation 

et + u- ve + A^e ^0, « = curi A^'^e, 

where V is the spatial gradient operator and A = (— A)^/^. The unknown function 
is 0: [0,T] X R2 -> R and we assume a G (0, 2]. 

2. Error bounds for Strang splitting: statement of results 

Let r be a positive integer and associate with it 

p = r + 2£-l, q = r + e-l^p-e, 

where £ > 2 is the degree of the polynomial P in (jl.ll) . The integers p, q and r will 
be kept fixed throughout the paper. 

We require local well-posedness of (|l.ip in H'~ and H'^ in the following sense: 
For both s = r and s = q and for a given time T, there exists R > such for all 
Mo G iJ" with ||uo||_f/= < ^7 there is a unique strong solution u G C{[0,T], H'') of 
(II. ip with initial data uq, and the dependence on the initial data is locally Lipschitz- 
continuous: There is a constant K — K{R,T) < oo such that the solutions u, u 
corresponding to arbitrary initial data uq, uo in the iJ^-ball of radius R satisfy 

(2.1) \\u{t)-u{t)\\„. <K\\un-ua\\H^ for < t < T. 

For the KdV equation, it is known from [T3] that this well-posedness assumption 
holds for every s G N, for arbitrary final time T and with arbitrary positive R. 

For the general class (|l.ip . (|1.2p . it can be shown with arguments as in Lemma [X^ 
below, that the above well-posedness result holds, at least, with sufficiently small 
R = R{T) for any given T, or for sufficiently small T = T{R) for any given R, 
the same as for the inviscid Burgers equation. For the equations listed in the 
introduction, however, the well-posedness situation is more favorable than for the 
inviscid Burgers equation; see [H] (the viscous Burgers equation) , ^ [T31 [H] (the 
KdV equation), [3] (the Benney-Lin equation), and [11] (the Kawahara equation). 

We consider solutions bounded by 

(2.2) \\u{t)\\Hp<p<R forO<i<r. 
In particular we assume that the initial data are in iJ^: 

uo e HP. 

Under these assumptions we will in this paper show the following results for the 
Strang splitting (fT^ . 

Theorem 2.1. (First-order convergence in H"^ and bound in HP) There is 
Ai > such that for At <Ai and tn = nAt < T, 

\\Un - u{tn)\\Hi < C'lAt and \\Un\\HP < Cq. 

Here, At, Cq and Ci only depend on \\uq\\hi', P, and T. 

This result will be used to prove our main result: 

Theorem 2.2. (Second-order convergence in H^) Assume that we have a 
solution u of (|l.ip that satisfies (j2.2p . Define the Strang approximation u„ hy 
(II. 8p . Then the following statement holds: There is At > such that for At < At 
and tn — nAt < T , 

\\un-u{tn)\\Hr < CaAt'. 
Here, At and C2 only depend on \\uq\\hp, P, and T . 
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For the KdV equation, this result yields second-order convergence in H"^ for 
initial data in H'^'^^ for arbitrary r e N. This improves significantly on the result 
of [8] , where second-order convergence is shown in the H^-woiva for r > 8 for initial 
data in H^^^. The proof uses a regularity result for the inviscid Burgers equation as 
in [5], but is otherwise quite different from the proof given there. It is conceptually 
closely related to the error analysis of Strang splitting given in |17J for the nonlinear 
Schrodinger equation. The proofs of Theorems 12.11 and 12 . 2 1 will cover the remainder 
of the paper. We remark that throughout our computations, C denotes a generic 
constant whose value may change at each occurrence. 

3. Regularity results for the inviscid Burgers equation 

The following variant of a result in [8 will play a key role in the proof of Theo- 
rem 

Lemma 3.1. // \\<^bM\\hi < a for <t < At, then ||$*3(uo)||hp < e'="*||uo||Hp 
for < i < At, where c is independent of uq and At. 

Proof We find that w{t) = <^bM satisfies (cf. [8]) 



II ^ II II ^ ^ U^t / \ l|2 



2 dt 



(3.1) 



p 



3=0 



dlw dl{wwx) dx 



1.1. 

j=Q k=0 

For j < p, any of the above terms can be estimated by 



{w,Wt)HP 

diwdt^^wdt^wdx. 



diwd'l^^wdi ^wdx 



<C\ 



gnia.x{k+l,j~k} ^ 



L2 



Qmin{k+l,j-k} ^ 



L2 



\HP 



< Ca , 



by using the Sobolev inequality ||w||i=o < ^ (so that C = 1/V2) and the 

fact that minjfc -t- 1, j — fc} < ■^{j + 1) < ■^P < P — £ ^ q since p > 2£. For j ^ p 
we estimate for k < q and fc g — 1 by 



/ := 



dPrndt+^wdP-'^wdx 



< WQMl^ \\d: 



,k+l 



W\ 



w\ 



L2 



<C\\w\\ 



\HP \\W\\Hk+2 \\uJiiBp-k ■ 

To estimate this further we have to distinguish the cases k + 2 < q and k ^ q. In 
the former case we find 

I < C\\w\\^p \\w\\h, \\w\\hp < Calkll^p , 
while in the latter case we obtain 

IHP II "^llHi - ^" II "'llz/p ' 

£ < £ + r - 1 = q. 

\- 1 ~ q or q + 2 < k + 1 < p we estimate 



w\ 



I<C\ 



w\ 



HP 



W\ 



since q + 2<q + £ = p and p — q - 
In the two cases where either k 



w\ 



:k+l 



w\\l2 \\dl~ 



<C\\w\\hp\\w\ 



Hp ll"'llif'« + i II ""llffp-fc + i • 

In both cases this is bounded by Ca ||i«||^p, since in the case fc + 1 = q we have 
p-k + l= £ + 2<2t<p, and for fc + 1 > + 2 we have p - k + 1 < p - q = £ < q. 
We are left with the term where k = p — j , viz. 

1 
2 



d^w d^.'^^ww dx 



{d^w) dxW dx 
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<C\\w\\^, \\w\\%^ , 



since q>£>2. Thus, 

d 
It 

which concludes the proof. □ 



\w{i)\\H^ < ca\\w[t)\\^p , 



We also need the following lemma: 

Lemma 3.2. // ||uo||ir™ < M for some m > 1, then there exists t(M) > such 
that ||$^(Mo)||ff'" < 2M for < t < t{M). 

Proof. The calculation (|3.ip with m instead of p, and the subsequent arguments, 
show that 

and hence the result follows by comparing with the majorizing differential equation 
y' = cy"^- □ 

Finally, we will make use of the following regularity result. 

Lemma 3.3. // ||ito||/fp < Cq, then there exists t depending only on Cq, such that 
the solution of the inviscid Burgers equation with initial data uq, w{t) — $^(wo), 
satisfies 

w € C2([0,t],i?«) and weC^{[0,t\,H''). 
Proof. If we define, for i G [0,t\ with t from Lemma 13.21 for m = p, 

w{t) ^uo+tB{uQ)+ [ {t~ s)dB{w{s))[B{w{s))]ds 
Jo 

and not4!l dB{w)[B{w)] — w'^Wxx + ^wwl, then we have w € C^([0, t], i/'). Since 
= dB{w)[B{w)] = B{w)t and t?t(0) = B{uo) = wt{Q) and w(0) = uq = w(0), 
we have w = w. The second statement is proved in the same way, by computing 

Wttt- □ 

4. Local error in H'^ 
Lemma 4.1. The local error of the Strang splitting is bounded in H'^ by 

\\^^\uo)-^T+bM\\h. <ciAt', 
where Ci only depends on \\uo\\hp- 

Proof. The proof follows that of the local error bound in [171 §4-4]. We write 
e^'^v = ) to indicate the linearity of the flow of A. We start from the variation- 
of-constants formula for u{t) = ^l^_^_g{uo), 

(4.1) u{t) = e''^uo+ f e^'-''^'^B{u{s))ds, 

Jo 

which is just the formula ^(i) — ^(0) — (p{s) ds for (f>{s) — e*^*^'*^'^M(s), and its 
sibling formula 

(4.2) B{u{3)) = B{e'^uo) + f dB[e'^'-'''^'^u{a))[e^'-'''^^ B{u{(j))] da, 



We adopt the notation from [T]. 



HOLDEN. LUBICH, AND RISEBRO 



which is nothing but the formula B{ip{s)) — B{ip{0)) — dB{ip{a))[ip{<7)] da for 
ip{a) = e(^"'^)^M(CT). We insert (g^) into (gJ]) with i = At to obtain 

u{At) = e^'^uo + / e(^*-")^B(e"^uo)ds + ei 
Jo 

with 

/•At i-s 

(4.3) ei=/ / e(^*-*)^dB(e(*-'^^^u((T))[e(^-")^B(w((7))]dCTds. 



On the other hand, for the result after one step of the Strang splitting, 

we use the first-order Taylor expansion with integral remainder term in i?*, 

(4.4) =v + AtB{v) + At^ f\l - 9)dB{^%^\v))[B{<^%^*{v))] d9. 

Jo 

This is justified for v — e^'^^/^uo <= and for sufficiently small At by Lemma 
We therefore obtain 



AtA 

ui — C 

with 



u, = e—uo + Ate'^'^/^B{e^'^^\o) + 62 



(4.5) 62 = At^ /V - 0)e^*^/2^B($ff *(e^*^/2yo))[S($|''*(e^*^/'Mo))] d0. 

Jo 

The error thus becomes 

pAt 

(4.6) ui~u{At)^Ate'^'^/^B{e^'^^^uo)- / e(^*-^)^B(e^^uo) + (e2 - ei), 

and hence the principal error term is just the quadrature error of the midpoint rule 
applied to the integral over [0, At] of the function 

(4.7) fis) = e(^*-^)^B(e^^wo). 

We express the quadrature error in first-order Peano form. 



pAt pi 

At/(iAt)- / /(s)ds = At^ / Ki 
Jo Jo 



(d) r{9At)dd, 



where ki is the real-valued, bounded Peano kernel of the midpoint rule. We find 

with the Lie commutator 

[A,B]{v) = dA{v)[Biv)] - dB{v)[Av] = P(9,)(wx) - {P{d,)v)v, - vP{d^)v,. 
We have that P is linear in for k < £, and we compute 

Hence the terms containing derivatives of order £ + 1 disappear, and we obtain the 
commutator bound 

\\[A,B]{v)\\h. < C{Ml^\M\hp + \\v^\\l^\\v\\hp + \\v\\ji,-i) < C\\v\\l,. 

Since e*'^ does not increase the Sobolev norms, it follows that 

\\f{s)\\H.<C\\u,\\%,, 
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and hence the quadrature error is 0(Ar) in the i/* norm for uq G H'\ The H"^ 
norm of the remainder term 62 — ei is bounded by CAt for uq £ for sufficiently 
small At (by using Lemma [5T^ . We include the details for the convenience of the 
reader: 



ei 



Jo Jo 

((e(--)^z.(a))(e(--)-45(z.(a))))^ 
|lB(u(cr))|l^,+i dads 



da ds 



At 



Jo 

At rs 



Hi 

da ds 



Jo 

At rs 



< c 

< c 

< CAt^R 



l'"(f^)llff<!+i ll'"(f^)llff<!+2 d(^ds 



Jo 

2 p3 



and 



< Af 



< At 



(\i-e) 

Jo 



,AtA/2 



de 



HI 



< At'C C $|^*(e^*'^/2^o) B($|^*(e^*^/2^o)) 

Jo 



de 



< AfC 



if 9+2 



^6* 



de 



< CAt^R^. 
Thus the proof is complete. 



5. Proof of Theorem 12.1 



□ 



The proof uses "Lady Windermere's fan" with error propagation by the exact 
flow [51 p. 160, Fig. 3.1]. In our approach the necessary regularity for estimating 
local errors by Lemma 14.11 is ensured by Lemma 13.11 via an induction argument, 
which we present next. 

We make the induction hypothesis that for fc < n — 1, 

\\uk\\m < R, \\uk\\Hp < e^^'"'^'\\uo\\Hp < Co 
\\uk ~ u(tk)\\Hi < "fAt, 



where Cq 



_ „2cRT\ 



uoll-f/p with c from Lemma 13.11 and 7 = K{R,T)ci(Cq)T 



with K{R,T) from the local Lipschitz bound (|2.ip and ci(Co) is the constant of 
Lemma [13] for starting values bounded by Co in H'p . We then show that the above 
bounds also hold for fc = n as long as nAt < T and At is sufflciently small. 
We denote, with <i>* = for brevity, 

which is the value at time i„ of the exact solution of p.ip starting with initial data 
Uk Sit time tk- We note 

Un = U^, U{tn) = U° . 

We estimate 

n-l 



\\Un-u(tn)\\Hi < ^ 



k+1 ku 



k=0 
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n-1 

fe=0 

For fc < n — 2 we have ||\I'^*(ufe)||^f<! — ||ufe+i||H<! < R, and 

||$^*K)||h. < ||<i>^*K)-$^*(u(ife))||H. + ||$^*(«(ife))||H. 

< K{R, At)\\uk - u{tk)\\H^ + ||u(ifc+i)||ff, 

< K{R,At)"fAt + p, 

which is bounded by R if At is so small that 

K{R,At)jAt < R- p. 
Using (|2.1|) and Lemma [4. II we therefore have, for fc < n — 1 and nAt < T, 

<X(i?,T)||^P^*K)-$At(y^)||^, 

< K{R,T)ci{Co)At\ 
With this estimate we obtain, again noting nAt < T, 

\\un-u{t^)\\H. <nK{R,T)ci{Co)At^ <jAt. 
With ^At < R ~ p, we have 

||Wn||H<! < -R- 

Since ||$^(-y)||HP < \\v\\hp, Lemmas O and O yield, for At <t{R), 

This completes the proof of Theorem 12.11 

6. Local error in 
Lemma 6.1. The local error of the Strang splitting is bounded in H'' by 

where ci only depends on \\uq\\hp- 

Proof. The proof follows the lines of [17j §5.2]. Instead of (|4.4p we now use the 
second-order Taylor expansion 

= " + AtB(w) + ^At^dB{v)[B{v)] 

+ dB{<^>%^\v))[dB{<^>%^\v)m'^>%^\v))]]) d0 
where henceforth we abbreviate the integral remainder term as 

At^ f {d^B{B, B) + dB dB B^ (^^'^'(w)) dO. 

Hence, 

m = e'^^^uo + Ate^*^/25(e^*^/2«o) + ^At^e^'''^^dB{e^*''^^uo)[B{e^'^/\o)] 
+ At"" f W~ efe^'"^'^ (d^B{B, B) + dB dB ($|^*(e'^*^/2^o)) dO 

where 62 is given by (|4.5p . 
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In (|4.3p we express the integrand by a formula of the type (|4.2[) by using 

G{u{a)) = G(e"^Mo) + T dG{e^'''^^^u{T))[e'^^-^^^ B{u{t))] dr 
Jo 

for the last part of the integrand in (|4.3p , with 

Giv) = G.^^iv) = di?(e(^-^)^z;)[e(^-'^)^i?(«)], 

and 

dG{v)[w] =d^B{eS'-''^^v) [e('*-'")^w,e(^-'^)^B(w)] 
Then we obtain, 



ei 



At 



/O ^0 "'0 

We return to the error formula (j4.6p and write the principal error term 

<.Ai 

Ai 6^*^/25(6^*^/2^0) - / 6(^*-^)^B(6^^wo)ds 

Jo 

in second-order Peano form 

Atfi^At) - / /(s) ds = At^ / KzCfi*) f'iOAt) d9 
Jo Jo 

with the second-order Peano kernel K2 of the midpoint rule and / of (|4.7p . We have 

The double commutator reads 

[A, [B,A]]{v) = [A,dAB{v) - dB{v)A{v)]{v) 

= dA^B(v) - dAdB{v)A{v) - d^AB{v)A{v) 

- dAdB{v)A{v) + d^B{A{v)f + dB{v))dAA{v) 
= A\B{v)) - 2A{dB{v)A{v)) + d^B{A{v)f + dB{v)A^{v) 
= A^{vv.,) - 2A{vAv), + {{Avf), + {vA^v),, 

since d'^A ~ 0. The operator A is linear in 9^, and the highest order derivative in 
the term {{Avf )^ \s t + 1 <2l-\, so we consider 



V. 



We note that in the above sum, the terms containing derivatives of order 2£ + 1 
and 2£ disappear, so 

\\f"mH^<G\\uo\\i.. 

Now for the difference of and (|i?5|) . 

62 — ei — ^At'^g{^At,^At) — / / g{s,a) da ds + 62 — ei, 

Jo Jo 

where 

g{s,a) = e(^*"^)^dB(6"^uo)e("-"'^B(6"^uo) 
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and 



and 



ei 



62 



At ps pa 



dG. 



Jo Jo 



drdads. 



At^ 5(1 - f )^e'^*'^/2 (d'^B{B, B) + dB dB b) (^^'^'(uo)) dO. 



To estimate the remainder terms e, we calculate 



dG,^^(w)u;||^. < Id's (e(^-")^B(w),w) ||^_ + ||dB (e'^^'^^w) e("^")'^dB(w) 
< C(||i?(t;)||^.+i ||u;||^„+i + \\v\\^,.^, \\dB{v)w\\H^^,) 



< C 



ff '-+2 



< C||i;||^,+2 ||w||ff.+2 , 



and 

\\{d'^B{B,B)+dBdBB) {v)\\^, < \\d'^B{B{v),B{v))\\^^ + \\dB{v)dB{v)B{v)\\ ^ 

< C (\\B{v)\\l.^, + ||t;||^„+, \\dBiv)B{v)\\ ^r- 



<C\\v\\%.+s. 

Then, using Lemma |3. 11 the remainder terms are bounded by 

PiIIh- + l|e2|lH. < CAtm\uo\\%, + hoIlD, 

since 2i-l>3. 

The first two terms in 62 — ei are the quadrature error of a first-order two- 
dimensional quadrature formula, which is bounded by 



iAt^5(iAi, iAt) - / / g{s,a)dads 
Jo Jo 



max ||9(7/ 9s||//r -I- max \\dg/d(j\\H'-) , 

where the maxima are taken over the triangle {(s,ct) : < ct < s < At}. In order 
to estimate the partial derivatives we write 

g{s,cr) = e^^*-''>'^dB{v{s))w{s,a), 

where 



v{s) = e^^Mo 



and w{s,a) = e^'-^^^Bivia)). 



With this notation 

^ = e(^t-s)A (^AdB(v(s))w(s, a)) + d"^ B(Av(s),w(s, a)) + dB(v(.s))Aw(s, a)) 

OS 

= e^"^*"*^^ {-A{v{s)w{s, <t)) + Av{s)w{s, a) + v{s)Aw{s, a))^ . 

Since the flow determined by A does not increase the Sobolev norms, it suffices to 
estimate the H'^'^^ norm of —A{vw) + {Av)w + vAw. We see that the derivatives 
of order £ vanish, so that 

Og 2 3 

< C ||w(s)||^,+j \\w{s, <T)\\jjr+e < C \\uo\\ Hr+i \\uo\\ H-+e+i < C \\uo\\ , 



ds 



since i>2. For the other partial derivative, we get 



da 



= e^^'-'^^dB{v{s)) (e^'-"^^ {-AB{v{s)) + dB{v{a))Av{(j))) , 
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SO that 

dg 



da 



<C\\v{s)\\h^^, \\-AB{v{s)) + dB{v{a))Av{a)\\^^+, 

The last factor is —A{vvx) + {vAv)x, and again the derivatives of order ^+1 vanish, 
so that 



dg 



da 



< C'lkoll^p ■ 



Therefore we obtain 

||e2-ei||^. <CAi3(||uo||HP + ||?^o|lD, 
which together with the bound for the quadrature error of the midpoint rule for / 
yields the stated result. □ 

7. Proof of Theorem 12.21 

We use "Lady Windermere's fan" in H"^ in combination with the bound of m„ 
from Theorem 12.11 and the local error bound from Lemma 16.11 For t„ = n/S.t < T 
we have, as in the proof of Theorem 12. 1[ 

n-l 

fc=0 
n-l 

< ^if(i?,r)||vi/^*K)-$^*K)|U. 

fc=0 
n-l 

< ^(^' r)c2(Co)Ai3 < K{R, T)c2{Co)TAt\ 

k=0 

This completes the proof of Theorem 12.21 
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